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ABSTRACT 

In this paper we present a parameter estimation analysis of the polarization and temperature power spectra 
from the second and third season of observations with the QUaD experiment. QUaD has for the first time 
detected multiple acoustic peaks in the E'-mode polarization spectrum with high significance. Although QUaD- 
only parameter constraints are not competitive with previous results for the standard 6-parameter ACDM cos- 
mology, they do allow meaningful polarization-only parameter analyses for the first time. 

In a standard 6-parameter ACDM analysis we find the QUaD TT power spectrum to be in good agreement 
with previous results. However, the QUaD polarization data shows some tension with ACDM. The origin of this 
l-2cr tension remains unclear, and may point to new physics, residual systematics or simple random chance. 
We also combine QUaD with the five-year WMAP data set and the SDSS Luminous Red Galaxies 4 th data 
release power spectrum, and extend our analysis to constrain individual isocurvature mode fractions, constrain- 
ing cold dark matter density, a c d m i < 0.11 (95% CL), neutrino density, a n di < 0.26 (95% CL), and neutrino 
velocity, a nv i < 0.23 (95% CL), modes. Our analysis sets a benchmark for future polarization experiments. 

Subject headings: CMB, anisotropy, polarization, cosmology 



1. INTRODUCTION 

The anisotropies of the Cosmic Microwave Background 
(CMB) radiation are among the most important tests of cos- 
mology. The large-angle Sachs-Wolfe effect, multiple acous- 
tic oscillations and the Silk damping tail in the tempera- 
ture power spectrum have now been confirmed by a range 
of experiments from the largest angul ar scales dow n to an - 
gular scales of a few arcminutes (Dunk lev et al.l (2008); 
iReichardt etail (l2008h ). The full repository of CMB data 
available, in conjunction with other cosmological observ- 
ables, such as data coming from the large-scale distribu- 
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tion of galaxies or Supernova type la observations, are ex- 
tremely well described by the spatially fiat ACDM cosmolog- 
ical model. 

A generic prediction of cosmology is that the CMB pho- 
tons should be polarized at the 10% level. The polarization 
field can be decomposed into two components: primary even- 
parity, curl-free E'-modes are generated at the last scattering 
surface by both scalar and tensor metric perturbations (grav- 
itational waves); primary odd-parity B-modes are generated 
only by tensor perturbations due to gravitational waves pass- 
ing through the primordial plasma. Secondary anisotropies 
in both the E and B-mode polarization arise at the epoch of 
reionization, while E and fi-modes are mixed by gravitational 
lensing by in tervening large-scal e structure along the line of 
sight (see eg iHu & W hite (19971)). Observations of this lin- 
early polarized component provide an important consistency 
check of the standard model and a detection of primordial 
gravitational waves in the odd-parity B-mode on large angular 
scales would be strong evidence f or inflation. 

After the DASI 1 experiment dKovac et al.l (l2002h ) made 
the first measurement of Zs-mode power, other experi- 
ments have provided us with furt her measurements a t 
a wide range of angular scales dBarkats et al.l (T2 005): 
iReadhead et al.l (12004]): iMontrov et all (120061): ISievers et al.1 
(12007 1 ) : iPage et all (120031) : iBischoff et all (120081) : iNolta et all 
(2008)). Despite this we were still lacking precision measure- 
ments of its power spectra down to arcminute scales, as we 
have for the temperature. The B-mode polarization has not yet 
been detected and only upper limits have been determined. 

The QUaD 2 experiment is at the forefront of this small- 
scale polarization quest, and after three years of observations 
has delivered the highest resolution /s-mode spectrum and the 

1 DASI stands for "Degree Angular Scale Interferometer". 

2 QUaD stands for "QUEST and DASI". In turn, QUEST is "Q & U 
Extragalactic Survey Telescope". The two experiments merged to become 
QUaD in 2003. 
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tightest upper limits on the B-modes yet measured. This is a 
significant improvement over the first se ason of data results, 
previously reported by lAde et all (120081) . In particular, the 
sensitivity of QUaD has allowed us to see, for the first time, 
four acoustic oscillations in the E-mode spectrum and all sig- 
nificant oscillations in the TE spectrum to I = 2000. The over- 
all consistency of peak phases and spacin gs between the tem- 
perature and QUaD EE data was shown in lPrvke et al.l (120081) 
(hereafter referred to as the "Power Spectra Paper"). 

In this paper we concentrate on using the QUaD tempera- 
ture and polarization power spectra to constrain the standard 
cosmological model. Using this baseline model, we analyze 
the different contributions coming from each of the QUaD 
spectra. We also go beyond the standard LCDM model us- 
ing the QUaD data in combination with WMAP and SDSS to 
constrain an isocurvature contribution. 

2. COSMOLOGICAL PARAMETER ESTIMATION: METHODOLOGY 
2.1. Monte Carlo Markov Chain 

The Monte Carlo Markov Chain (MCMC) method is a 
method designed to efficiently explore an unknown Probabil- 
ity Distribution Function (PDF) by sequentially drawing sam- 
ples from it according to a proposa l probability function in 
our case the Metropolis algorithm (Metrop olis et alJ (119531) 
among others). The ensemble of these samples constitute a 
Markov Chain whose distribution corresponds to that of the 
unknown PDF. We adopted the Gelman and Rubin /^-statistic 
to verify that our chains are properly mixed and converged 
dGelman & Rubin! dl992l) ; IVerde et al.1 (120031) ). and compare 
sets of 4 chains of around 100 000 steps from which we re- 
move a burn-in period during which our criterium is not met. 
We use at least 80 000 steps after burn-in. 

The rate of convergence of the Markov chain is slowed 
down by degeneracies between parameters, and the choice of 
the step size in the Metropolis algorithm. Therefore we apply 
a standard p artial re-parametrization of the parameter space as 
suggested in Kosow sky et al.1 ([2002 ). To further reduce the re- 
maining degeneracies between paramet ers, we apply a change 
of basis in parameter space described in Teg mark et al.l (2004) 
which uses a covariance matrix to take account of all the cor- 
relations between parameters. In the new basis, the new pa- 
rameters have zero average and unit variance. 

When we have a fair sample of the underlying distribution, 
the MCMC method trivializes marginalization to a simple 
projection of the points of the chain. The mean marginalized 
value of each parameter will hereafter be called "the mean re- 
covered model". Note that some authors refer to the mean 
recovered model as the "best-fit model". 

To obtain constraints on the mean parameters values one 
simply produces the ID histograms of the chain values for 
each parameter, and calculate confidence interval s using the 
p' h an d (l-p th ) quantiles of the histograms as in lVerde et alJ 
( 120031) . Normally we use 68% equivalent to a nominal l-cr 
constraint. However if the constraint as defined above hits the 
prior boundary on one end we instead choose the level which 
contains 95% of the total probability and quote an upper (or 
lower) limit. 

We shall also plot 2-D marginalized parameter distributions 
with 68% and 95% contours estimated at AlnL = -2.3 and 
—6.17 from the peak values. Assuming a Gaussian distribu- 
tion we quote \ 2 values corresponding to our mean recov- 
ered model, and the Probability To Exceed (PTE), P(> ^\v), 
which gives the random probability to have found the mea- 
sured value of x 2 or greater by chance, for v degrees of free- 



dom. 

2.2. The Likelihood and Nuisance Marginalization 

The likelihood for the measured Q bandpowers is well ap- 
proximated by a Gaussian distribution, given by 



P(Q,|0<xexp 



(1) 



where ACb(p) = Cb~ C b h (p), Cb are the measured QUaD 
binned Q bandpowers, C l b (p) are the theoretical power spec- 
tra which depend on the cosmological parameters, p, and 
which have been transformed to predictions of the binned 
spectra by means of the experimental bandpower window 
function (BPWF), andM w = (AQACy) is the measured C b 
bandpower covariance matrix (BPCM). 

The BPCM is estimated from an ensemble of simulations of 
the CMB sky, assuming a fixed fidu cial ACDM, run th rough 
the QUaD analysis pipeline (see e.g. iBrown et alJ (120051) ). In 
our analysis, this BPCM remains independent of the cosmo- 
logical parameters, but in principle we should vary it as we 
move around parameter space. This re-scaling of the sample 
variance with C th is equivalent to saying that P(C l b (p) | C b ) 
is log-normally distributed, and in the c ase where noise is 
present, offset log-normally distributed. iBond et all d!998l) 
propose a transformation of the bandpowers and BPCM 
which accounts for this effect and allows use of a fixed 
BPCM. We have tested the impact of this transformation on 
our full data set and find that the difference in parameter esti- 
mates and uncertainties is insignificant. 

The Gaussian likelihood has the added benefit that we can 
simplify the marginalization over beam and calibration nui- 
sance parameters. Indeed, in this case, o ne can directly appl y 
an analytic marginalization scheme as in Brid le et al.1 ((2002). 
The resulting marginalization results in extra terms added to 
the BPCM of our Gaussian likelihood that acts as a source of 
extra noise. The likelihood function is then given by 



InL = -- ACtM'bl AC£ - ^TrlnM', 



where 



M' bb , = M bv + a^CbCl + 2a 2 b 5£ 2 b Cb£ 2 b,Ct, , 



(2) 



(3) 



is the marginalized bandpower covariance matrix, where a^ai 
is the variance on the calibration, a b = 0FWHM/voln2, Qfwhm 
is the effective beam size, S is the fractional beam error and 
lb is the average multipole in a bin. 

In assessing the goodness-of-fit of our mean recovered 
models, and when comparing the measured Cb bandpowers 
with the WMAP5 ACDM model, we shall use the x 2 -statistic, 
introduced previously in Section lXTl and which we define by 



X 



'-ACbM'llAcl, 



(4) 



where we use the nuisance marginalized BPCM, M' b b' , as de- 
fined in the previous equation. 

2.3. The Standard Cosmological Model 

We parameterize our flat ACDM cosmological model with 
the following standard set of 6 cosmological parameters: the 
Hubble constant, Hq = lOOh kms _1 Mpc _1 ; the physical matter 
density, Sl m h 2 \ the physical baryon density, O^/z ; the ampli- 
tude of scalar fluctuations, A s ; the scalar spectral index, n s ; 
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QUaD Temperature and Polarization 
QUaD Polarization 
QUaD Temperature 
WMAP 



FIG. 1. — 2-D projected basic parameter likelihood surfaces with two-parameter 1- and 2-sigma contours for QUaD only constraints using the TTITEIEEIBB 
data set (TP: the blue contours), using the TEIEEIBB data set (P: the red & magenta contours) and using the TT spectrum (T: the yellow & orange contours) 
versus the WMAP5 constraints (black/empty contours). Pivot scale used is k„ = 0.05 Mpc -1 . 



and the optical depth, r. When using QUaD data by itself, we 
present the combination A s e~ 2r as our individual constraints 
on the degenerate parameters A s and r are prior driven and 
thus biased, as explained in Appendix|A] Initial conditions are 
taken to be purely adiabatic with an initial power-law mass- 
density perturbation spectrum. Due to the range of angular 
scales probed by QUaD, the pivot point we use when ana- 
lyzing QUaD data by itself is k p = 0.05 Mpc -1 (note this is 
independent of h). When comparing our QUaD results with 
WMAP we regenerate WMAP best-fit values based on this 
pivot value using our own pipeline, however when adding 
QUaD data to other data sets for a combined analysis we re- 
vert to the WMAP preferred pivot scale of k p = 0.002 Mpc" 1 . 

To generate our the oretical spectra w e use the publicly 
available CAMB code dLewis et al.1 d200Qt) X including the ef- 
fects of reionization, and gravitational lensing by foreground 
structure. We impose the following flat priors in the likeli- 
hood analysis: < Q, c h 2 < 1, < Vt h h 2 < 1, 0.005 < 9 < 0.1, 
< r < 0.8, < A s < 2.5 and < n s < 2. The parameter 9 is 
the angular sound horizon, and fl c h 2 is the physical cold dark 
matter density. Note that the pa rtial re-parametr i zation of the 
parameter space as suggested in lKosowskv et al.l (|2002) intro- 
duces an implicit prior on the h parameter. 

3. RESULTS: BASIC 5-PARAMETER CONSTRAINTS 

3.1. QUaD Only Constraints 

The QUaD data set we use to constrain cosmological pa- 
rameters are the optimally combined spectra obtained from 
the 100 GHz, 150 GHz and frequency-cross temperature and 
polarization E- and B- power spectra, measured in 23 band- 
powers over angular multipoles from 200 < £ < 2000 as de- 
scribed in the Power Spectra Paper. We estimate a 10% un- 
certainty on both the calibration (in power units) and the beam 
sizes, assuming an effective beam of 4. 1 arcminutes 3 . 

3 The data set is publicly available online at http://quad.uchicago.edu/quad 



In parameter estimation we use the diagonal and the first 
two off-diagonal terms of the BPCM for TT-TT, TE-TE, EE- 
EE, and BB-BB covariances, but only the diagonal and first 
off-diagonal terms in the TT-TE and TE-EE covariances. This 
is motivated by the need to avoid excessive noise in the off- 
diagonal terms of the BPCM, due to its estimation from nu- 
merical simulations. We also ignore the covariance between 
TT and EE, which is much smaller than the other terms. 

We explore various sets of combinations of the QUaD tem- 
perature and polarization data in order to understand the new 
information each spectrum brings to parameter estimation. In 
Table Q] we present the mean recovered models. Figure Q] 
shows the corresponding 2-D marginalized contour projec- 
tions of the likelihood in the 5-parameter space. 

All of the statistics shown verify that the QUaD TT tem- 
perature power spectrum is compatible with the results from 
WMAP5. This is a non-trivial test, since the overlap of scales 
measured by QUaD and WMAP5 is only in the range £ ss 200 
to £ w 950, while the QUaD data extends to £ w 2000 with 
good signal-to-noise. 

However the analysis involving the polarization spectra 
are in less good agreement, yielding a high value of the 
baryon content. Indeed for the TEIEEIBB combination we 
have n b h 2 = 0.03 19 ±0.0046, compared to Q b h 2 = 0.02261 ± 
0.00062 from WMAP5. The X 2 of the WMAP5 best-fit model 
for this data set has a PTE of 7% indicating a modest degree 
of tension. The spectrum responsible for this behaviour seems 
to be the TE spectrum. 

Clearly the TE only constraints are weak and most param- 
eters are prior driven, but surprisingly we obtain constraints 
on £li,h 2 and VL„,h 2 that are not influenced by their choice of 
priors. Figure [2] shows the 2-D projected likelihood surface 
for the (ilbh 2 ,il m h 2 ) parameter space. To illustrate the dif- 
ference with the TT only contours we overplot them. In ad- 
dition we show the results from WMAP5 and the Big Bang 
Nucl eosynthesis (BBN) co nstraint of Q. b h 2 = 0.0214 ±0.002 
from lKirkman et all (12.003b . We also show in Figure|3]the TE 
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TABLE 1 

Basic cosmological mean parameter constraints using QUaD bandpower spectra for various data 

combinations. 



Symbol Q08 TT/TE/EE/BB Q08 TEIEEIBB Q08 TT Q08 TE Q08 EEIBB WMAP5 



n b h 2 0.0334 0.0319 ±0.0046 0.0242 0.0398 ± 0.0051 0.0366 ™j>!^ 0.02261 ^™^ 

n,„h 2 0.119™'™ 0.117±0.017 0.133 iopff 0.149±0.025 0.155 5{J|f 0.1329 ^Kg 

h 0.91 ±0.09 0.90±0.10 °- 75 5'l8 0.87±0.11 0.77±0.18 0.717 

A s e- lT a 0.66±0.08 0.63±0.09 0.64 +° : ^ 0.63±0.13 0.79±0.21 0.614 ™ : °j| 

«/ 0.809±0.078 0.766±0.152 0.848™|'J 1.337 +°g' 0.534^{^ 0.967±0.015 



\ 2 {uf 88.60(86) 74.78(63) 12.73(17) 19.67(17) 33.16(40) 

PTE: P(>x 2 W) 40.26% 14.72% 75.38% 29.14% 76.94% 

X 2 (WMAP5|208) C 108.63(92) 86.99(69) 14.48(23) 31.44(23) 41.62(46) 

PTE(WMAP5\Q08) 11.36% 7.07% 91.24% 11.24% 65.60% 

a The pivot point for A s and n s is k p = 0.05 Mpc" 1 for both the QUaD data and WMAP5 data. 

b x 2 f° r the 6-parameter mean recovered model against QUaD data, with the number of degrees of freedom in brackets. 
c x 2 f° r WMAP5 mean recovered model given the QUaD data set, with the number of degrees of freedom in brackets. 



0.25 




0.00 0.02 0.04 0.06 

n b h 2 

FIG. 2. — 2-D marginalized contours of the parameters f2(,/i 2 versus 
£l m h 2 obtained from QUaD TE data only. Also plotted are the contours 
from QUaD 7Tdata o nly, the results from WMAP5, and the BBN constraint 
tKkkman et aT] <200l ). 

QUaD bandpower spectrum versus its mean recovered model 
and the WMAP5 full data set best-fit model. This figure vi- 
sually illustrates the differences of the mean recovered mod- 
els between the two data sets in terms of height and location 
of the peaks. The main reason for the higher baryon density 
parameter seems to be larger acoustic oscillations at higher 
multipole in TE, as well as a shift to higher multipoles of the 
peaks, resulting in a slight degeneracy with h, which may ex- 
plain its high value. The origin of this source of tension is 
unclear, but could be due to a new physical mechanism, resid- 
ual systematics or random chance. 

Another interesting result comes from the EE and BB spec- 
tra. As expected they provide very little information on pa- 
rameters, and the x 2 of the WMAP5 best-fit model for this 
data set has an acceptable PTE of 65.6%. They do however 
have an unusual feature; the preferred range of scalar spectral 
index values is low (n s = 0.534 ±o i6i)- S° although TT and TE 
share the majority of the constraining power, the EE and BB 
spectra exert an influence in combination with the remaining 
spectra by restricting the «,-range to low values. 

If we combine all the spectra together then the polarization 
data dominate the constraints. The majority of parameters are 
consistent with the WMAP5 results, but the spectral index n s 
is lower, influenced by the EEIBB contribution, and the Hub- 
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FIG. 3. — Plot of the QUaD TE data bandpower spectrum (in red) versus 
the QUaD TE mean recovered model (in solid black) together with WMAP5 
best-fit model (blue). For the TE mean recovered model, we assumed the 
WMAP5 best-fit value for the optical depth (r = 0.087), which corresponds 
to an Amplitude of A s = 0.75, given our A s e~ 2l ~ constraint. 

ble parameter h and ttbh 2 are higher, driven by the polariza- 
tion data in particular the TE spectrum. Compared to the BBN 
value ilbh 2 is almost 3-cr away. Figure [4] shows a comparison 
between our TT/TE/EE/BB data and models and the WMAP5 
data and best-fit model. The \ 2 of the WMAP5 best-fit model 
indicates that there is an 11.36% chance that the combined 
QUaD spectra are a realization of this model. 

3.2. Combining QUaD with Other Data Sets 

In this section we will add to the QUaD spectra the WMAP 
5-year TT, TE, EE and BB data set. We use the WMAP5 like- 
lihood code, publicly available on the LAM BDA website 4 , 
and their methodology (Dunkley et al. (2008)), but do not in- 
clude the Sunyaev Zel'dovich (SZ) marginalization. We shall 
assume the two data sets are independent, as the QUaD data 
only covers a small fraction of the WMAP5 sky, and the over- 
lap in multipole range is only partial. 

We will further add large-scale structure data from the 
SDSS Luminous Red Galaxies (LRG) fourth data release 
using publicly avail able likelihood code , measurements and 
window functions (Tegm ark et al.l (l2006l) \ Results from the 
SDSS LRG and the main SDSS galaxy samples are consis- 
tent, but the former provides higher signal-to-noise ratio. We 

4 LAMBDA website: http://lambda.gsfc.nasa.gov/ 
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TABLE 2 



10000 
8000 
6000 
4000 
2000 



Basic mean parameters for QUaD TTITEIEEIBB, SDSS LRG and WMAP5 data. 


Symbol 


Q08 TT/TE/EE/BB+WMAP5 


Q08 7T/r£/££/BB+WMAP5+SDSS 


WMAP5 


n b h 2 
n m h 2 

h 

T 


n n?9^ +0.00055 

U.UiZjj _ .00057 

0.1266 ±0.0060 

0.733 ±0.027 

0.087 ±0.017 

0.805 ±0.038 

960 +0 014 
u.you _ 013 


02235 -rt- 00051 
u.u^jj -0.00058 

1266^™™ 
U.1ZOO .0 0039 

0.731 ±0.019 
0.087 ±0.016 

0.960 


02259 + 00061 
_ oo62 

0.1329 ±0.0065 

71 i +0.027 
u - ' 1J -0.026 

0.087 ±0.017 
0.816±0.039 
966 +0 014 



1 The pivot point for A., and n s is k p = 0.002 MpcT 1 for QUaD, WMAP5 and SDSS LRG data. 





















QUoD 2008 - 










WMAP 5yr _ 










WMAP CIs _ 














500 


1000 


1500 


2000 2500 




1000 1500 
multipole, I 




1000 1500 
multipole, I 



' ■ ■ » i i a i 



500 1000 1500 2000 2500 

multipole, 

FIG. 4. — The QUaD power spectra used in this analysis (red points) for 
TT, TE, EE and BB (top to bottom). The blue data points are WMAP 5- 
year power spectra data. The blue line shows the basic WMAP5 best-fit 
model as defined in Dunklev et al. 12008), while the black solid line shows 
our TTITEIEEIBB mean recovered model with values given in Table [T] We 
assumed the WMAP5 best-fit value for the optical depth (r = 0.087), which 
corresponds to an Amplitude of A , = 0.79, given our A s e r constraint. 

use the SDSS LRG matter power spectrum over wavenum- 
bers 0.07 /iMpc" 1 < k S Dss < 0.2/iMpc" 1 so that we do not have 
to consider any nonlinear correction. We marginalize over 
the amplitude of the galaxy power spectrum which removes 
any dependence on the galaxy bias parameter, b g , and linear 
redshift-space distortion. 

As can be seen in Table [2] the QUaD TTITEIEEIBB power 
spectra have little impact on the baseline 6-parameter mean 
parameter fit when combined with WMAP5. This is per- 
haps unsurprising, given the accuracy of the WMAP5 mea- 



surement of the first acoustic peak in the TT spectrum, and 
its low-£ power in TT and TE. The impact it does have is to 
tighten the error bars on parameters determined from the rela- 
tive heights of acoustic peaks, i.e. on the baryon density, fi^/z 2 , 
and the matter density, Vt m h 2 , as QUaD data adds a substantial 
amount of well-defined peak information at high-f. 

When we combine the SDSS LRG and WMAP5 data with 
the QUaD data we see an improvement compared to the 
QUaD and WMAP5 combination, as expected. This improve- 
ment is mostly due to the extra constraining power on £l m h 2 
and VLbh 2 coming from the break-scale in the SDSS galaxy 
power spectrum, and the baryon acoustic oscillations. How- 
ever the QUaD data still reduces the uncertainty on flbh 2 . 

4. BEYOND THE STANDARD 6-PARAMETER MODEL: 
ISOCURVATURE MODES 

Theoretical predictions of isocurvature modes and their 
evolution, and the role of CMB polarization observa- 
tions in constraining them , has been an acti v e fiel d 
over the past f ew ye ar s ([K awasaki & Sekiguchi (2007); 
Keskit alo etail (l2007[) ; iBean et all (l2006l) : iBeltran et all 
(2004); Moodle v et alj (l2004l) among many others). Pure 
isocur v ature perturbations hav e been ru led out (IStompor et alJ 
dl996l); lLanglois & Riazuelol (l2000l) ; lEnqvistet alf ^OOO); 
Amendola et al. (2002)) although the presenc e of a subdom- 
inant i socurvature fraction has been claimed ( Keskital o et alJ 
(12007b '). Observationally, isocurvature modes have a phase 
difference from adiabatic modes, which provides a distinct 
signature. 

We can completely characterize the primordial perturba- 
tions by one adiabatic and several isocurvature components. 
The adiabatic component is given by the associated curva- 
ture perturbation 7Z corresponding to an initial overdensity 
S = 8p/p. The non-adiabatic components are given by en- 
tropy perturbations S x = <5 A -(3/4)(5 7 between photons and a 
different species, x. These correspond to four possible nonde- 
caying isocurvature modes: baryon density, cold dark matter 
densit y (cdmi), neutrino density (ndi) and neutrino velocity 
(nvi). iBucher et al] d2000) have presented a thorough analysis 
of these components. 

We parameterize the contribution of adiabatic and isocur- 
vature modes to the total temperature and polarization power 
spectra by 



Cf =A 2 [(l-a)tf M +aCf' u °] 



(5) 



where a is the isocurvature fraction. The adiabatic spectra, 
C^' Ad , and the isocurvature spectra, C^' Iso , are defined with 
unit amplitude and the same spectral index. In this analysis 
we shall assume there is no correlation between adiabatic and 
isocurvature modes, and will constrain one isocurvature mode 
at a time. Also we do not present results for the baryon density 
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TABLE 3 

CDM ISOCURVATURE MEAN PARAMETER CONSTRAINTS FOR QUAD TTITEIEEIBB, WMAP5 AND SDSS LRG DATA. 



Symbol 


Q08 TT/TE/EE/BB+WMAP5 




Q08 TT/TE/EE/BB+WMAP5+SDSS 


WMAP5 


Q b h 2 


fl (IMP +0.00080 
u.uzjiz n.00081 




0.02280 ± 0.00070 


02362 + 00( » 6 
u.uzjoz a nno94 

1279 ^ 

0.759 ±0.041 


Cl m h 2 
1, 


a ]2 i 4 +6.0067 
u,u " -0.0068 

0.773 ±0.039 




1256 +° 0038 

u ' lzJO -0.0037 

0.746 ±0.022 


T 


0.087 ±0.017 




0.084 


0.087 ±0.017 


A/ 

a cdml (95 % cl) 


0.786 ±0.037 
0.987 ±0.023 
< 0.19 




0.797 |g 

a g 76 +0.018 
u y/ ° -0.017 

< 0.11 


0.789 ±0.038 

A AAO +0.027 

<0.21 


a The pivot point for A., and n s for all isocurvature constraints is k p 


= 0.002 Mpc" 1 for QUaD, WMAP and SDSS data. 






TABLE 4 




NDI Isocurvature mean parameter constraints for QUaD TTITEIEEIBB, WMAP5 and SDSS LRG data. 


Symbol 


Q08 TT/TE/EE/BB+WMAP5 




Q08 77777i/££/fiB+WMAP5+SDSS 


WMAP5 


tt b h 2 


0.02370 ±0.001 10 




0.02300 ± 0.00080 


02410 +0 00130 

1270 +00075 
u.iz/u _ 0077 

780 +°- 058 
U. /OU _o o5 5 

0.090 ±0.017 


VLji 1 
h 

T 


1200 + 00073 
u.izuu _q 0072 

0.800 ±0.030 

0.090 ±0.017 




0.1261 ±0.0039 

a 7 ci +0.027 
0.085 ±0.016 


A s 


0.869 ±0.058 




A oc<; +0.053 

-0.051 

976 +™ 
-0.017 

<0.26 


0.872 ±0.059 


n s 

a nil (95 % cl) 


a one +0.0271 
•j.yyj _o.()252 
<0.37 




1 003 +0031 

i.uuj _ 0030 

<0.38 



TABLE 5 

NVI Isocurvature mean parameter constraints for QUaD TTITEIEEIBB, WMAP5 and SDSS LRG data. 



Symbol 


Q08 TT/TE/EE/BB+WMAP5 


Q08 TT/TE/EE/BB+WMAP5+SDSS 


WMAP5 


n b h 2 
n m h 2 


0.02350 ±0.00090 
0.1260 ±0.0060 


07339 +0 00080 
u.Uijj7 -0.00070 

0.1277 ±0.0040 


0.02390 ±0.00100 
0.1 330 ±0.0063 


h 

T 


0.745 ± 0.029 
0.088 ±0.017 


0.734 ±0.018 
0.087 «jg« 
0.854 


a 790 +0.027 
U - //!S -0.028 

0.089 ±0.018 


A s 


°- 851 S 


0.859 ±0.048 


n s 

a nvl (95%cl) 


0.980 »;i 
<0.27 


0-978 S 
<0.23 


0.988 
<0.27 



isocurvature mode as these have a very similar signature to the 
cold dark matter mode. 

We analyze the QUaD TTITEIEEIBB power spectra com- 
bined with WMAP5, and combined with the WMAP5 plus 
the SDSS LRG data. The shape of the galaxy power spec- 
trum is sensitive to an isocurvature contribution, and has 
been used in the past to impro v e on isocurvature constraints 
(e.g. lBeltran etail (l20Ml2005h : ICrottv et all (l2003h l. There- 
suits we obtain are given in Table [3] for the cdmi mode, in 
Table|4]for the ndi mode and in Table[5]for the nvi mode. 

Our analysis shows a small improvement in the isocurva- 
ture cold dark matter constraint when we add the QUaD to 
the WMAP5 data, from a cdmi < 0.21 to a cdmi < 0.19 (95% 
confidence limits). In addition we find an improvement in the 
tlbh 2 and il m h 2 constraints. There is a similar improvement 
for the neutrino density isocurvature constraints: we go from 
«ndi < 0.38 to a n & < 0.37. For the neutrino velocity isocur- 
vature modes there is no improvement, the constraint staying 
at a nv j < 0.27. 

We can further improve on these results by adding the SDSS 
LRG data. When we do this the cold dark matter isocurva- 
ture constraint becomes a c d m i < 0.1 1 (95 % CL). The largest 
improvement is for the neutrino density isocurvature mode, 



«ndi < 0.26, while the smallest improvement is for the neu- 
trino velocity isocurvature mode, a nv \ < 0.23. 

5. SUMMARY AND CONCLUSIONS 

We have presented a standard cosmological parameter con- 
straint analysis, and its extension to include isocurvature 
modes, using QUaD 7T, TE, EE and BB bandpower spectra. 
This is the first CMB experiment to detect with confidence the 
acoustic oscillations in the EE spectrum and so the first to be 
able to provide significant constraints on cosmological param- 
eters from individual CMB polarization spectra — in fact our 
polarization-only constraints are superior to our temperature- 
only constraints. In combination with the WMAP5 data set 
QUaD offers a small improvement in the constraints on the 
baryon and matter densities. 

We find our QUaD temperature data is in good agreement 
with the results from WMAP5, which is a non trivial test 
of LCDM as the QUaD data extends to I w 2000 with good 
signal-to-noise. However, our polarization (TE, EE and BB) 
data is in less good agreement yielding a higher baryon den- 
sity value of n b h 2 = 0.03 19 ± 0.0046, compared with VL b h 2 = 
0.0242 from our TT data and 0.02261 from our 

re-analysis of WMAP5. A \ 2 test shows there is a 7% proba- 
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FIG. 5. — 2D marginalized plot, showing the scattered values (green) of the 5 parameters mean recovered b asic cosmological model obtained from 50 
simulations of QUaD TT, TE, EE and BB data generated from the WMAP3 best-fit model (in red) from table 2 of Spergel et al. 12007). The mean over the 50 
simulations is shown as a blue point. 
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bility of the QUaD polarization results arising by chance, as- 
suming the WMAP5 ACDM model is correct. Although not 
of high significance, this modest level of tension, that seems to 
originate from the TE spectrum, could be due to new physics 
in polarization, residual systematics in the data, or random 
chance. It will be interesting to see if this trend continues in 
future polarization experiments. 

We also investigate isocurvature cold dark matter den- 
sity, neutrino density and neutrino velocity modes. We find 
QUaD provides a marginal improvement on the fractional 
cold dark matter density mode parameter, a C dmi, from < 0.21 
for WMAP5 alone to < 0.19. 

QUaD is funded by the National Science Foundation in 
the USA, through grants AST-0096778, ANT-0338138, ANT- 
0338335 & ANT-0338238, by the UK Science and Technol- 



APPENDIX 

SIMULATING PARAMETER ESTIMATION 

We test our MCMC pipeline by running it on a set of 50 simulations of the QUaD TT, TE, EE and BB bandpower spectra. 
These were generated by simulating the signal and noise properties of the time-ordered data and passing these through the QUaD 
pipeline in the same way as the data (see the Power Spectr a Paper for de t ails). T he input cosmological model for these simulations 
was the WMAP3 mean recovered model (see Table 2 of Sper gel et alj ([2007)). The scatter in the values of the mean recovered 
model obtained from each one of the 50 simulations can be seen in Figure [5] We also overplot the average values for each 
parameter calculated from the 50 simulations (see blue points). We have verified that the scatter in the simulated mean parameter 
results is close to the size of the contours produced by our MCMC code when using real data, indicating that our code accurately 
estimates the parameter uncertainties. 

We can also compare the mean parameter values and the scatter about them with the input WMAP3 best-fit model (red crosses). 
The average over the simulations closely matches the input model indicating that our parameters are not biased. If constrained 
independently, the scalar amplitude, A s , and the optical depth, r, parameters are biased, their values being systematically higher 
than the input values. This is due to the combination of the large degeneracy between the amplitude, A s , and optical depth, r, and 
the parameter priors. To break this degeneracy requires large-scale polarization measurements probing the re-ionization bumps 
at lower £-modes. As can be seen in the Figure, this problem can be avoided if we combine A s and r into the parameter A s e 



along the line of degeneracy, which is the approach followed in Section [3T 
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